DEPOLARIZATION  DIFFUSION  DURING  WEAK  SUPRATHRESHOLD 
STIMULATION  OF  CARDIAC  TISSUE 

Vladimir  Nikolski,  Aleksandre  Sambelashvili,  and  Igor  R.  Efimov, 
Department  of  Biomedical  Engineering,  Case  Western  Reserve  University,  Cleveland, 

Ohio  44106-7207,  USA 
E-mail:  ire@cwru.edu 

Abstract 

Background:  The  theory  of  stimulation  of  cardiac  tissue  based  on  virtual  electrode  (VE) 
polarization  was  presented  by  Wikswo  and  his  group  (Wikswo,  Biophys.  J.,  1995).  Four 
basic  regimes  were  distinguished:  anodal  make,  anodal  break,  cathodal  make,  and 
cathodal  break.  We  sought  to  detennine  applicability  of  this  theory  for  near-threshold 
stimulus  intensities. 

Methods:  For  the  experiments  we  used  fluorescence  optical  imaging  to  map  virtual 
electrode  polarization  at  the  epicardium  of  Langendorff-perfused  rabbit  hearts  stained 
with  voltage-sensitive  dyes.  Unipolar  cathodal  pacing  was  applied  during  diastole  via  a 
platinum-iridium  electrodes. 

Results:  During  near-threshold  cathodal  make  stimulation  we  found  deviations  from 
Wikswo's  theory.  We  observed  that  after  quick  virtual  electrode  polarization  activation 
starts  at  hyperpolarized  regions  and  then  causes  a  subsequent  activation  of  the 
depolarized  regions.  Such  an  activation  pattern  appears  similar  to  break  activation.  The 
effect  of  the  depolarization  diffusion  from  depolarized  to  hyperpolarized  regions  of  the 
VE  obtained  by  us  experimentally  was  confirmed  theoretically  by  2D  active  bidomain 
simulations. 

Conclusions:  For  weak  near-threshold  stimuli  the  activation  pattern  does  not  fit 
completely  into  the  established  model  of  Roth  and  Wikswo  thus  revealing  its  limitations. 

Introduction 


Ability  of  electric  point  stimulation  to  produce  a  response  in  excitable 
tissues  [1,2]  and  induce  [3]  or  terminate  [4]  arrhythmia  in  the  heart  is  well  known.  Yet, 
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the  exact  mechanisms  of  electric  stimulation  have  been  unknown  until  a  recent  discovery 
of  virtual  electrode  polarization  produced  by  point  stimulation  [5,9]  which  results  in  a 
peculiar  "dog-bone"  pattern  of  positive  and  negative  polarizations.  These  polarizations  of 
opposite  signs  are  the  so-called  virtual  cathode  (VC)  and  virtual  anode  (VA).  They 
represent  the  driving  force,  which  can  be  mathematically  expressed  as  the  activating 
function  [10,11].  The  last  is  governed  by  two  major  parameters:  the  gradient  of 
extracellular  electric  field  and  the  structural  heterogeneity  of  the  heart,  both  contributing 
to  polarization  of  the  cellular  membrane  during  stimulus.  Active  ionic  properties  of  the 
heart,  particularly  Ca-channels,  modulate  these  polarizations  [12]. 

Dekker  [13]  demonstrated  that  both  the  onset  (make)  and  termination  (break)  of 
stimulation  of  appropriate  intensity  and  duration  could  produce  a  propagated  response. 
Roth  [14]  and  Wikswo  et  al.  [9]  provided  the  first  mechanistic  explanation  of  the  "make" 
and  "break"  excitation  (ME  and  BE)  based  on  virtual  electrode  polarization.  They  showed 
that  during  make  stimulation  the  heartbeat  originates  at  virtual  cathode(s)  and  is  delayed 
at  virtual  anode(s).  In  contrast,  during  BE  depolarized  regions  are  unexcitable  due  to 
inactivation  of  sodium  channels  by  the  long  stimulus.  Therefore,  during  break  stimulation 
the  heartbeat  originates  at  virtual  anode(s)  and  is  delayed  at  virtual  cathode(s). 

The  goal  of  our  study  was  to  investigate  the  mechanisms  of  stimulation  with 
clinically  relevant  short  stimuli  of  weak  near-threshold  intensities. 

Methods 


Experimental  methods 

The  experiments  were  performed  on  perfused  cardiac  preparations  obtained  from 
New  Zealand  rabbits  (n  =  12).  The  heart  was  placed  onto  a  Langendorff  apparatus,  where 
it  was  retrogradely  perfused  with  oxygenated  modified  Tyrode's  solution  as  previously 
described  [15]. 

The  hearts  were  stained  by  0.5  pM  of  voltage-sensitive  dye  di-4-ANEPPS. 
Motion  artifacts  in  optical  recordings  were  suppressed  by  excitation-contraction 
uncoupler  2,3-butanedione  monoxime  (BDM)  (15  mM). 


The  heart  was  paced  at  a  steady  state  cycle  length  of  400  ms  at  the  base  of  left 
ventricle  with  Basic  Pacing  Electrode.  Test  stimuli  were  applied  with  Test  Stimulus 
Electrode  in  the  middle  of  field  of  view  at  a  coupling  interval  of  350  ms  with  amplitude 
of  0.2-60  mA  and  duration  of  2  ms  from  a  constant  current  source  (A385,  WPI). 
Recording  Electrode  was  used  to  detect  the  propagated  response  in  order  to  confirm 
suprathreshold  pacing.  Basic  pacing  stimulus  strength  was  adjusted  to  twice  the  diastolic 
threshold  of  excitation. 

Transmembrane  potential  (mV)  was  recorded  using  a  16x1 6-element  photodiode 
array  (Hamamatsu)  with  high  resolution  (252  pm,  343  ps)  from  the  4x4  mm  area  near  the 
tip  of  the  test  electrode  (0.12  mm  platinum-iridium  Teflon-coated  wire)  positioned  on  the 
epicardium  of  the  left  ventricle. 

Numerical  methods 

We  guided  our  experiments  with  theoretical  predictions.  We  reproduced  2- 
dimensional  bidomain  model  with  active  kinetics  as  described  by  Trayanova  and 
Scoubine  [18]. 

The  bidomain  model  is  based  on  the  representation  of  the  cardiac  tissue  as  two 
interpenetrating  extra-  and  intracellular  domains  each  of  them  having  different 
conductivities  along  and  across  the  direction  of  the  fibers  [16,17].  The  state  variables 
describing  the  system  are  intracellular  (j)i  and  extracellular  (f)e  potentials  defined 
everywhere  in  the  domain  of  interest  Q .  The  variable  of  interest  is  the  transmembrane 
potential  defined  as  a  difference  Vm  —  (f)i  —  (f)(, .  The  following  coupled  reaction-diffusion 
equations  constitute  the  bidomain  model: 

v-(tt  Mi)  =  im 

V )  =  -/,„ -/0  inQ 

where  and  (7  are  intra-  and  extracellular  conductivity  tensors  respectively,  I m  is  the 
volume  density  of  the  transmembrane  current,  70  -  volume  density  of  the  stimulation 
(shock)  current. 

The  transmembrane  current  is  described  as  a  sum  of  capacitive,  ionic  and 
electroporation  currents  [18]: 


4  =  P(Cm  -zf  +  444,0  +  <44,0  ■  4) 

ot 

where  [3  -  the  surface  to  volume  ratio  (total  membrane  area  divided  by  the  total  tissue 

volume),  Cm  -  specific  membrane  capacitance,  G(Vm,t )  -  electroporation  conductance. 

The  last  is  described  by  empirical  equations  [18],  e.g.: 

dG  ,  o 

—  =  aexp(/?(4 -VrestY){\ -exp(-r(4 -Vrest)2)),  G(0)  =  G0 

dt 

Here  OC,/3,y  are  electroporative  coefficients,  Vrest  is  the  resting  transmembrane 

potential.  The  values  of  all  the  parameters  are  given  in  Table  1. 

As  a  model  of  cardiac  myocytes  we  used  the  original  BRDR  model  [19]  modified 
by  Skoubine,  Trayanova  et.  al.  in  order  to  accommodate  strong  electric  fields  [18]. 

We  used  the  following  boundary  conditions 

M  =  0  (j)  =  0  on  <3Q 

on 

For  numerical  solution  of  the  bidomain  system  we  used  a  square  grid  with  the 
space-step  hx=hy=0.4  mm,  the  time  step  was  At=0.02  ms.  To  invert  the  sparse  matrix 
inversion  we  employed  the  Generalized  Minimum  Residual  (GMRES)  method  [20,21]  as 
the  most  robust  and  fast  for  our  case.  We  used  diagonal  preconditioning  for  this  method. 
Calculations  were  performed  on  Dell  Pentium  III  PC. 

Results 

During  excitation  by  strong  anodal  and  cathodal  stimuli  at  the  first  5  ms  the 
anodal  wavefront  has  an  elliptical  shape  elongated  along  the  fibers,  while  the  cathodal 
wavefront  is  more  rectangular.  The  activation  starts  at  the  positive  VE  regions  with  a 
delayed  propagation  through  hyperpolarized  area.  These  results  correspond  to  those 
previously  obtained  by  Wikswo  et  al.  for  so-called  make  stimulation. 

However,  the  near-threshold  current  strength  stimulation  develops  differently. 
The  anodal  wavefront  is  formed  in  a  similar  way  to  the  one  for  the  strong  stimulus  while 
cathodal  activation  starts  from  initially  hyperpolarized  regions  but  not  from  the  positive 


VE  areas  (Fig.  1).  As  a  result  cathodal  wavefront  for  near  threshold  stimulus  originates 
and  propagates  alike  an  anodal  wavefront. 

The  observed  difference  in  a  wavefront  fonnation  during  strong  and  near 
threshold  cathodal  stimulations  remains  in  the  absence  of  excitation-contraction 
uncoupler  2,3-butanedione  monoxime.  Without  the  uncoupler  the  pacing  threshold  is 
significantly  smaller  so  it  is  difficult  to  the  see  the  VE  pattern  for  the  near  threshold 
stimulus.  Nevertheless  the  comparison  of  the  initial  wavefronts  indicates  that  during  the 
near  threshold  cathodal  stimulation  wavefront  also  starts  from  the  hyperpolarized  area  of 
VE. 

Computer  simulations  confirmed  our  experimental  observations.  Near-threshold 
2  ms  stimulus  applied  to  the  center  of  the  square  sheet  of  bidomain  tissue  resulted  in 
virtual  electrodes  and  after  termination  of  the  stimulus  the  wavefront  started  at  the 
initially  hyperpolarized  regions.  This  is  illustrated  in  Fig  2  where  two  action  potentials  at 
two  different  points  of  VE  are  presented.  As  it  can  be  seen  the  virtual  cathode  area 
depolarizes  first.  Yet  the  depolarization  is  not  strong  enough  for  the  activation  and  it 
diffuses  away  to  the  hyperpolarized  regions  of  the  virtual  anode.  Those  regions  fire  in 
response  to  the  intensive  current  inflow  and  then  their  activation  spreads  back  to  the 
virtual  cathodes  causing  them  to  fire. 


Conclusions 

This  work  reveals  certain  limitations  to  the  commonly  accepted  model  of  make 
and  break  excitation  suggested  by  Wikswo  et  al.  Particularly  it  shows  that  for  weak 
stimuli  depolarization  diffusion  leads  to  a  typical  break-like  excitation  pattern  for  the  case 
when  make  excitation  is  expected. 

The  nature  of  the  ionic  mechanisms  that  underlie  the  interplay  between  excitation 
and  depolarization  diffusion  and  their  role  in  wavefront  formation  still  remains  obscure. 
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Parameter 

Value 

1-0  /uF / cm2 

p 

3000  cm~l 

^ix 

0.375  mS/cm 

°iv 

3.750  mS/cm 

aex 

2.140  mS/cm 

^  ey 

3.750  mS/cm 

a 

2.5  •  10  3  mS/cm2  -  ms 

P 

2.5  - 10  5  \/mV2 

Y 

1.0  - 10  9  \/mV2 

h 

50  fiA/cm 2 

Table  1.  2D  active  bidomain  model  parameters 


Fig  1.  Right  panel  -  experimental  recording  of  VE  at  2  ms  after  the  near  the  threshold 
make  cathodal  pulse.  Left  panel  -  action  potential  traces  from  hyper-  (solid)  and 
depolarized  (dashed)  areas.  The  hyperpolarized  regions  activate  before  the  depolarized 
regions. 


ms 


Fig  2.  Right  panel  -picture  of  VE  at  2.4  ms  after  the  near  the  threshold  make  cathodal 
pulse  obtained  from  the  active  bidomain  simulation.  Left  panel  -  action  potential  traces 
from  hyper-  (solid)  and  depolarized  (dashed)  areas.  The  hyperpolarized  regions  activate 
before  the  depolarized  regions. 


